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We study the effects of periodic matter in kaon regeneration, motivated by the possi- 
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^\j ■ bihty of parametric resonance in neutrino oscillations. The large imaginary parts of 

t^ ! the forward kaon-nucleon scattering amplitudes and the decay width difference AF 

prevent a sizable enhancement of the Kl — > Ks transition probability. However, some 
interesting effects can be produced using regenerators made of alternating layers of 
^ , two different materials. Despite the fact that the regenerator has a fixed length one 

CL|| can obtain different values for the probability distribution of the K^ decay into a final 

o ■ state. Using a two-arm regenerator set up it is possible to measure the imaginary 

D . parts of the i^^(^'')-nucleon scattering amplitudes in the correlated decays of the 

i?i)-resonance. Combining the data of the single-arm regenerator experiments with di- 
rect and reverse orders of the matter layers in the regenerator one can independently 
^ . measure the CP violating parameter 5. 
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Recently, there has been a renewed interest |jT[ in the possibihty of parametric resonance 
in neutrino oscillations in matter suggested in 0, ^. For a neutrino beam propagating in a 
medium with periodic density, one can obtain a large probability for the transition from one 
flavour state to another, even if the neutrino mixing angles both in vacuum and in matter 
are small. 

In nature, there are other systems similar to oscillating neutrinos, in particular the 
neutral mesons K^ — K^. Hence, it is interesting to investigate if one can obtain the 
parametric resonance in this case. The analogue of the neutrino weak flavour basis are K^ 
and K^ and the mass eigenstates are Kl and Ks- Since the former states are maximally 
mixed, it is obvious that one cannot enhance the K^ — K^ transition probability. However, 
in this case this is not the relevant question. 

Let us assume that we have a neutral kaon beam propagating in vacuum. After a time t 
larger than the Kg lifetime (r^ = 0.894 x 10~^° s) the beam is essentially a Kl beam. If this 
beam traverses a thin slab of material, a small Ks component will emerge, because K'^ and 
K^ have different scattering amplitudes. This is the well-known regeneration phenomenon 
(see, e.g., 0]). Assuming that the beam enters the regenerator at t = and denoting by 
\K^{t)) the state of the beam when it emerges, we have 

\Kii{t)) = {ks\T\KL)\Ks) + {kL\T\KL)\KL), (1) 

where {Ks,l\T\Kl) are the transition amplitudes in the regenerator, and {Ks,l\ are the 
reciprocal states (see Eqs. (||) and (^ below). If the regenerator is a medium with a density 
that is a periodic function of the coordinate along the beam direction, we would like to see 
if it is possible to enhance the Kl — > Ks transition amplitude. Our aim in this letter is to 
address this question. 

Assuming CPT conservation, but not CP conservation, the rest-frame evolution equation 
for the K^ — K^ system propagating in a medium is 

4 (f:ufa^?^>U !:;:). p) 



dt\K« J \ 1^ ^, + v j \K 

where t is the proper time (we follow closely the notation of ref. [Q). Hence, in vacuum 
{V = V = 0) the eigenstates of the Hamiltonian Hq are 

\Kl) = p\K') + q\K') (3) 

and 

\Ks) = p\K')-q\K^), (4) 

with the corresponding eigenvalues fiL,s = ^tT'L,s — -^^l,s, l^ = ifJ-L+fJ^s)/'^ and A/i = hl — jJ'S- 
Since the phase of p/q is of no physical significance, we shall assume this ratio to be realQ. 

"'^ Since the relative phases between p/q and certain ratios of ampHtudes do have physical content, the 
phase convention we use here implies a phase convention for the decay amplitudes as well. 



We write 

where (5 ~ 3 x 10^'^ ||^ is a measure of CP violation. Since CP is not conserved, the 
diagonahzation of Hq cannot be accomphshed with a unitary transformation. This, in turn, 
imphes the use of the reciprocal basis B 



(,?,, . i(i<A-vi(A-»|) 



(6) 



In a medium with Na nuclei per unit volume, V {V) is given in terms of the forward 
scattering amplitude /(O) (/(O)) for a K^ (K^) beam 0], i.e. 

V = -Mi/(0) , (8) 

m 

with the average kaon mass m = 2.52 fm~^. The simplest way to introduce a periodic 
medium is to consider two different elements with number densities Na and Nb positioned 
one after the other and to build a regenerator with k layers of this ab junction. The beam 
evolution through this multilayer regenerator can be described in terms of the evolution 
operator 

U, = UbUaX ... xUbUa, (9) 

"> V ■' 

K 

with 

Ua = exp{-tHata) (10) 

and 

Ub = exp{-iHbtb) , (11) 

where Hi {i = a, b) are the Hamiltoneans for layers a or b, given in Eq. (H). Since this is a 
2x2 matrix it is convenient to represent it using the Pauli a matrices. One can then write 

Ha = Fa + (TEa, (12) 

where 

Fa=^^+^^^, (13) 
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and Eq is a three dimensional vector with components 



Ei'^ - 


A^ /p ^ q\ 
4 \q p) ' 


E^^^ - 


,A/i fp q\ 
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Va-Va _ AK 
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which are complex numbers. Introducing the complex unit vector 

E„ 



(14) 
(15) 
(16) 



n, = -^, Ea = ^/Ea ■ Ea , (17) 

and 

^a = Eata , (18) 

one immediately obtains 

Ua = exp {-iFata) (cOS ipa - 1(7 ■ Ua siu ipa) , (19) 



With the obvious replacements a ^^ b one obtains from Eq. (p!9|) Ub- Then the product 
UbUa is 



UbUa = exp [-1 {FJa + Fbh)] [r - 20- ■ X] , (20) 



with 



Y = cos ipa cos ipb- sin ipa sin ipb{na- rib) , (21) 

X = sin (pa cos (pbiia + sin ipb cos (farib - sin (yj^ sin (pb{i^a x rif,) . (22) 

The vectors that we have introduced have complex components. However, the dot products, 
such as iIq ■ rifc, must be simply understood as 

n.-n, = X:n«n«. (23) 

i 

Notice that the third component of n^ x n^, is identically zero. Then Eq. (^) shows that 
X^^^ is symmetric with respect to the interchange of a and b. On the other hand, Eq. ([T5| ) 
shows that ni and n\^ vanish in the limit of CP conservation. Then, in this limit, X^"^^ is 
antisymmetric in a and b. Furthermore, in the same approximation the first component of 
Ua X lift is also zero. Hence X'-^^ is symmetric with respect to the interchange of a and b. 



A straightforward calculation shows that y^+X-X = 1. Then, defining another complex 
angle $ such that 



cos$ 
sin$ 



Y 



/X2 



it is possible to rewrite Eq. ( [2D| ) in the form 

Ui,Ua = exp [-i {Fata + Ffotf,)] exp ( -io--X $ 
with 



X 



X 



/X2 



(24) 
(25) 

(26) 

(27) 



This evolution operator is written in the K^—K^-hasis. Denoting Uab = UbUa, the symmetry 
properties of X*^*) deduced above enables us to obtain 



{K'\Uab\K') = {K'\Uta\K') 



but 



{K'\Uab\K^) - {K'\Uta\K'') oc S, 

i.e. the difference vanishes if CP is conserved. 

Finally, the evolution matrix for the propagation through k a6-layers is simply 



U^ = exp [-iK{Fata + Fbtb)] (cos(/t$) - icr ■ Xsin(K$) 



Inserting [/^ between the appropriate bra- and ket- vectors given by Eqs. (0)-(|5) and 
respectively, one obtains the Kl -^ Kl and the Kl — > Ks transition amplitudes 



and 



{Kl\U^\Kl) = exp [-m (F,t, + Fbh)] x 



cos(k$) sin(/T;$) 



(KsIU^IKl) = exp l-iK {Fata + Fbtb)] X 



X - sin(K<l>) { -2X3 + 
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(32) 
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(33) 



(34) 



Let us start our discussion with a careful examination of Eqs. (|14D-(|T^). The vectors 
Eq and Eft have the first components proportional to A/x/2 and the third components 
proportional to AV^/2. These quantities A/x and AV play a crucial role in the effect that 
we are searching for. On the contrary, the mean values n and [Vi + Vi)/2 are far less 
important. Their real parts disappear when we take the modulus square of the amplitude 
to obtain the transition probabilities and their imaginary parts give the overall damping 
factors. 

As a first approximation, we neglect CP violation. Let us further assume that A/x and 
AVi are real. A real A/x means that AF = 0. Although this is not true for the /C-meson 
system, there is no fundamental reason why it could not be so. Indeed, such a situation 
occurs closely for the B^ — B^ mesons. A real AVi implies equal imaginary parts for the K^ 
and K^ forward scattering amplitudes. As it is well known, this is not the case. This is in 
contrast with the case of neutrinos, where the absorption is weak, and to the leading order 
in weak interaction the scattering amplitudes are real. 

Within this unrealistic approximation it is possible to achieve a parametric resonance in 
Kl <-> Ks transitions in matter. The parametric resonance condition is X3 = ; we shall 
consider a particular realization of this condition in which the times t^ and tf, are chosen 
such that cos ipa = cos ip^ = 0. Then it follows from Eq. (|2^) that 



X = ±(n, X n,) . (35) 

As described above, the third component of the cross product n^ x n^ is identically zero 
and if one neglects CP violation the first component is also zero. In this approximation 
Eqs. (|33|) and ( ^4|) become 



(KlIU^IKl) = exp [-m (F,t, + Ffet,)] cos(fi:<l') (36) 

and 

(KsIU^IKl) = -exp[-m(F,4 + F,tb)]sin(fi:<l>)X2. (37) 

For an appropriate number of layers, k, one can suppress the Kl — ^ Kl probability and, at 
the same time, enhance the Kl — > Ks transition probability. To illustrate this effect, we 
plot the Kl -^ Ks transition probability as a function of k, in Fig. |l|. The calculation was 
done for a regenerator made of ^^Al and ^^''W and for an initial Kl beam obtained from the 
decay of the resonance at rest. The values of the K'^ and K^ scattering amplitudes on 
protons and neutrons were taken from ref . |^ . As we have explained, AF and the imaginary 
parts of AVi were set equal to zero. 

The times ta and tb were chosen in such a way that a complete Kl — > Ks conversion could 
be obtained. If we move away from this resonance condition we still obtain an oscillatory 
Kl —>■ Ks transition probability P{Kl —>■ Ks) but with a smaller maximal conversion. For 
instance, decreasing both ta and tf, by 17 % reduces the maximum value of P{Kl —>■ Ks) 
from 1 to 0.145. 



The resonance values of ta and tb (59.210 x 10~^^ s and 57.289 x 10~^^ s) are a factor of 
seven or six larger than ts- This by itself is sufficient to explain that the effect disappears 
as soon as we introduce the right values of F and AF, even with Ini(V^) = Ini(V^) = 0. We 
have checked that, in this case, P{Kl —>■ Ks) — 10~^ for n = 1 and decreases slowly with 
K. In addition, if we introduce the correct values for the imaginary parts of the scattering 
amplitudes, P{Kl — > Ks) for k = 1 is further reduced to 8 x 10~^ and even P{Kl —* K^) 
becomes 0.05, while in the previous case it was 0.93. 

Clearly, any measurable effect with kaons propagating in matter requires times of the 
order of r^. Unfortunately, for such times, even the toy model without imaginary parts 
gives a maximum value for P{Kl -^ Ks) of the order of 0.02 only. Then, the damping due 
to the imaginary parts washes out the effect. This is shown in Fig. ^ where we compare for 
the same ta and t^ P{Kl -^ Ks) in the toy model (Pi) and for real kaons traversing a real 
^'^Al - -"^^^W regenerator. 

So far, in all cases that we have considered, the total time that the particles spend in the 
regenerator, t = K{ta+tb), increases linearly with k. Obviously, after a few layers most of the 
particles will disappear due to their decay or absorption. Hence, it is interesting to examine 
another type of experiment, where the total time t is kept fixed, i.e. as k increases the times 
ta and tb are proportionally reduced. In Fig. |^ we plot P{Ki -^ Ks) as a function of k 
for this situation. For the beam velocity that we are considering, 10^^^ s corresponds to a 
pathlength of the order of 1 mm in vacuum. Then from Fig. ^ one can see that a regenerator 
made of a 12 mm layer of ^^Al followed by another layer of 12 mm of ^^^W (/t = 1) is less 
efficient than another regenerator with four alternating ^^Al - ^®^W layers of 6 mm each 
(k = 2). Perhaps this effect is better illustrated if, instead of the transition probability, we 
consider the decay of the kaons into a final state / after traversing the regenerator. From 
Eq. (|1]) one can calculate the time distribution P{Kfi(t) — > /) of the final state / after the 
kaon state initially produced as K^ passes through the regenerator and then spends outside 
it the proper time t (which for simplicity we took equal to the proper time spent inside the 
regenerator). The result is (e.g. ref. @]) 

P(KR{t)^f) = \{f\T\Ks)\'\{ks\U^m)\'x 

X [e-^st + |^^|2g-rit^2|wj|e-^*cos(%-Am-t)] , (38) 

where 

I I iSf _ {KL\Uf,\KL) 

{j^s\Uk\j^l) 

In our example, shown in Fig. |^, we have assumed that one measures the tt+tt" final state. 
The magnitude and the phase of ?7+_ were taken from ref. [^. The probability distribution 
after passing n (^^Al-^^^W) layer junctions increases with k. In the same Fig. ^ we also plot 
the probability distribution for a regenerator where the layers are in reverse order. In this 
case P{Kr -^ f) decreases with n, and both curves tend to a common limit. This is easy 
to understand. As the number of layers increase we are effectively approaching a "mixed 



material" with a density that has the average density of aluminum and tungsten. Since 
the regeneration effect is proportional to the density of the regenerator, one can understand 
that P{Kr — >■ TT^TT") increases with k for the ^^Al-^'^^'W regenerator and decreases in the 
^^^W-^^Al case 0. The variation with the order of the layers (notice that their total number 
is fixed) is a nice example of quantum mechanics interference. In this problem, the evolution 
matrix for each individual layer (cf. Eqs. ([T0| ) - ([TT1) ) is an element of the U(2) group. Hence, 
the evolution for the total number of layers is, of course, an element of U{2). Since U{2) is 
a non-Abelian group, shuffling the layers one obtains a different evolution operator. From 
this point of view. Fig. ^ is a consequence of the non-commutativity of the U{2) group. 

One should realize that the results shown in Fig. ^ are independent of the CP-violating 
parameter 5. However, it is possible to use this type of regenerators to measure CP violation 
at the 0-factories. To see how the effect arises let us recall that the 0-meson decays into the 
antisymmetric combination 

i= [\K\^))\K\-^)) - \K\p))\K\-p))] , (40) 

where p denotes the momentum of the particle. We assume that in the direction of — p 
we have a detector, called "left", and in the direction of p another detector, called "right". 
Both detectors measure muons from the semileptonic decays of the kaons. These decay 
amphtudes are 

(7r-/i+z/^|T|K0) = A, (41) 

{'K+^-v^\T\K^) = A*. (42) 

The kaons propagating to the right from the decay point have to traverse a regenerator made 
of two layers of different materials a and h. On the other hand, the kaons that propagate 
to the left must traverse a similar regenerator with two layers of the same width but in 
reverse order, h followed by a. With this setting one can show that the amplitude to detect 
in coincidence two /x"*" on both detectors is 

+ X [{K'\U,a\K'){K'\Ua,\K') - {K'^\U,a\K'^){K'\Ua,\K')]) . (43) 

We have introduced the ratio 

^ (vr>^.,|T|K°) 

between the AS" = AQ violating amplitude and the dominant one. Experimentally, x = 
[—2 ± 6 + i{1.2 ± 1.9)] X 10"^ 0, which is consistent with zero; theoretically, within the 

^For K = 1 the regeneration in the ^^Al-^^^'W case is less efficient than in the ^^^W-^^Al one because a 
fraction of Kl decays in aluminum before they reach a more efficient regenerator - tungsten. 



standard model one expects x ~ 10~^. Therefore in Eq. (|43|) we have neglected the term of 
order x^. Finally, let us point out that Eq. (|i3| ) goes trivially to zero when a=b. This is 
a simple consequence of the antisymmetry of the initial state. With a similar notation one 
can obtain the amplitude for two n~ in coincidence. The result is 

A{n\n^) = ^-^{{K'\Ua,\K'){K'\Ua,-U,a\K') + 

+ X* [{K'\Uba\K'){K'\Ua,\K') - {K'\U,a\K'){K'\Uab\K')]) . (45) 

The amplitude for a /i"*" in the left detector and a fi~ in the right detector is 

A A* 
A(7r-,7r+) = ^^{{K'\U,a\K'){K'\Uat\K')-{K'\Uta\K'){K'\Uat\K') + 
v2 

+ X* {K^\Uab\K^){K^\Uha-Uab\K^)+x{K'\Uab\K'){K^\Uab-Uta\K'')(^6) 

The other asymmetric amplitude A{7i^,7t~) is 

A A* - - 

A{n+,7r-) = -^{{K'\Uba\K'){K'\Uab\K')-{K'\Uba\K'){K'\Uab\K') + 

+ X* {K^\U,b\K'^){K^\Uba-Uab\K'>) + x{K'\UjK'){K^\Uab-Uba\K^)p) 

We shall now consider two asymmetries which can be measured in the two-arm experiments, 

R, ^ '^(-''-"^'^-l^^f'-^)'^ (48) 

|A(7r-,7r+)|^ 

' |A(7r-,vr+)|V|A(7r+,7r-)r 

and also an asymmetry which can be measured in the single arm experiments of the 
CPLEAR type (see e.g. ref. 0), 

K>(.-)i;-|-4U^^)i;. (50) 

Here 

A^bi^-) = A [{K^\Uab\K') + x{K'\U,b\K')] , (51) 

Abain^) = A* [{K'\Uba\K'') + x*(ir°|f/b„|K°)] . (52) 

The ratios Ri and R^ are CP-asymmetric observables. They depend on the intrinsic (i.e 
fundamental) CP violation parameter 6. Furthermore, since the regenerators are made 
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out of matter and not of equal amounts of matter and antimatter, they are themselves CP- 
asymmetric and so induce a macroscopic, extrinsic CP violation which in general contributes 
to both CP-violating observables, Ri and R3. However, interchanging the order of the layers 
leads to a partial cancellation of the extrinsic CP violating effects. For this reason the ratio 
i?3 is primarily sensitive to the fundamental CP violation. In the limit x = the cancellation 
of the extrinsic CP violation in R^ is exact and this leads to the result R^ ^ 26. This is 
not so for Ri which does not vanish when 6 = 0. The ratio Ri, for example, is normally of 
the order of unity as the extrinsic CP violation is of this order. For an aluminum-tungsten 
regenerator, using x = 0, ta = 24 x 10~^^ s and t^ = 12 x 10~^^ s, we find -Ri = 1.349 for 6 
given in ref. |^, whereas for 6 = the corresponding value is Ri = 1.334. One should notice 
that Ri is very sensitive to the imaginary parts of the effective Hamiltonian. Switching off 
the imaginary parts of the matter-induced potentials Vi and Vi reduces Ri by about a factor 
of 200, while switching off the decay rates F^ and F^, Ri would reduce it by about a factor 
of 8. If all the imaginary parts are set equal to zero, Ri is suppressed by a factor 2 x 10~^. 

The parameter R2 may appear as a CP-violating observable too, but in fact it is not. 
To see that one has to notice that CP transformation not only interchanges particles with 
their antiparticles but also flips the sign of all the momenta; for the two-arm setup under 
discussion this implies an additional interchange of the arguments of A{7ri,7rj) so that R2 
is unchanged under the CP transformation 6 -^ —6. It has a moderate sensitivity to the 
imaginary parts of the effective Hamiltonian. For the same regenerator and x = 0, we find 
i?2 = —0.69 with normal values of all the imaginary parts. Switching off Vi and Vi reduces 
|i?2| by about a factor of 2, while switching off F5 and T^ would reduce it by about a factor 
of 1.4. If all the imaginary parts are set equal to zero, R2 goes to zero. 

Thus, by measuring Ri and R2 one can obtain an information on the imaginary parts of 
the effective Hamiltonian of the K^K^ system in matter, and in particular on the imaginary 
parts of the K^ {K^)-imcleon scattering amplitudes. 

In conclusion, we have studied the effects of periodic matter in kaon regeneration. Mo- 
tivated by the possibility of the parametric resonance in neutrino oscillations in matter we 
considered similar effects in Kl — > Ks transitions. Unfortunately, the large AF and imagi- 
nary parts of the forward kaon-nucleon scattering amplitudes prevent a sizable enhancement 
of the Kl —>■ Ks transition probability (cf. Fig. |^). However, some interesting effects can 
be produced using regenerators made of alternating layers of two different materials. De- 
spite the fact that the regenerator has a fixed length one can obtain different values for the 
probability distribution of the Kl decay into a final state (cf. Fig. ^ and Fig. |^). Finally, 
we have pointed out that using a two-arm regenerator set up it is possible to measure the 
imaginary parts of the i^°(A'°)-nucleon scattering amplitudes in the correlated decays of 
the 0-resonance. Combining the data of the single-arm regenerator experiments with direct 
and reverse orders of the matter layers in the regenerator one can independently measure 
the CP violating parameter 6. 
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Figure 1: Kl — * Ks and Kl —* Kl transition probabilites P{LS) and P{LL), respectively, 
with lm{V) = lm{V) = AT = T = 0. The regenerator is ^^Al -^^^W and ta = 59.21 x 10"^^ s 
and U = 57.29 x 10~^^ s. 
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Figure 2: P2 is the K^ -^ Ks transition probabihty for ^^Al — ^^■^W regenerator with ta 
tf, = 8 X 10~^^ s. Pi is the same probabihty when neglecting all imaginary parts. 
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Figure 3: K^ — ^ Ks transition probability for ^^Al — ^^^W regenerator with t^ = t^ 
^ X 10-11 s. 




Figure 4: P{Kr -^ -n^ii )/(10^s ^) for two regenerators (Al-W or W-Al) with ta = h 
^ X 10-11 s. 
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